integer(kind=int4), parameter :: nLp = 15, np = 4
